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' Abstract. Due to the kinematic reversibility of Stokes flow, a body executing a 

, reciprocal motion (a motion in which the sequence of body configurations remains 

identical under time reversal) cannot propel itself in a viscous fluid in the limit of 
negligible inertia; this result is known as Purcell's scallop theorem. In this limit, the 
' \ Reynolds numbers based on the fluid inertia and on the body inertia are all zero. 

Q ' Previous studies characterized the breakdown of the scallop theorem with fluid inertia. 

O . In this paper we show that, even in the absence of fluid inertia, certain dense bodies 

undergoing reciprocal motion are able to swim. Using Lorentz's reciprocal theorem, we 
first derive the general differential equations that govern the locomotion kinematics of 
^ , a dense swimmer. We demonstrate that no reciprocal swimming is possible if the body 

motion consists only of tangential surface deformation (squirming) . We then apply our 
. general formulation to compute the locomotion of four simple swimmers, each with a 

I different spatial asymmetry, that perform normal surface deformations. We show that 

' the resulting swimming speeds (or rotation rates) scale as the first power of a properly 

defined "swimmer Reynolds number" , demonstrating thereby a continuous breakdown 
of the scallop theorem with body inertia. 
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1. Introduction 

At the small scales of swimming microorganisms, the inertial mechanisms of locomotion 
used by larger animals become ineffective. Fluid flow in the limit of zero Reynolds 
number is governed by the Stokes equations, which are both linear and independent of 
time. In this limit, Newton's equations of body motion reduce to balances of forces and 
torques, which depend on time only through the sequence of body configurations. Due to 
these properties, a periodic motion in which the sequence of body configurations remains 
identical under time reversal (termed a "reciprocal" motion) yields a net translation 
equal to its opposite, and therefore equal to zero. This result is known as Purcell's 
"scallop theorem" In particular, any periodic motion of an organism with one 
degree of freedom, such as a scallop, is reciprocal. Such an organism cannot propel 
itself in the limit of zero Reynolds number. 

Consider a swimmer (i.e., a self-propelled, deforming body) of density pp and 
characteristic size a immersed in a fluid of density p <^ pp and dynamic viscosity p. The 
swimmer executes a reciprocal motion of amplitude A and radian frequency u. There 
are three relevant Reynolds numbers that characterize this oscillatory motion The 
"unsteady" Reynolds number, Re^j = pa^u/ fi, is the scale ratio of the unsteady terms to 
the viscous terms in the Navier-Stokes equations. The "advective" (or "frequency" 0]) 
Reynolds number, Rcf = paAu/fi, is the scale ratio of the nonlinear advective terms to 
the viscous terms. Finally, the "particle" Reynolds number (or "Stokes number" [sl), 
Rcp = ppO^uj/fj,, is the scale ratio of the particle inertia to the viscous forces on the 
particle. 

How much inertial force is necessary for a reciprocal motion to become propulsive? 
Childress and Dudley j3] addressed this question for the case of a symmetric, reciprocal 
flapper and concluded that the breakdown of the scallop theorem occurs above a finite 
threshold of the "advective" Reynolds number, Rcf, of order unity. As discussed by 
Childress and Dudley, the breakdown of the scallop theorem is relevant to the study 
of swimming organisms that cross the Rcf threshold as theygrow, such as the larvae 
of certain molluscs 0], crustaceans [^], tunicats [gI, or fish [J 0]. The existence of a 
finite threshold above which the scallop theorem becomes invalid has been confirmed 
in laboratory experiments [ol, [13] and numerical simulations 11|. In contrast, Lauga 



[2I devised examples of oscillatory reciprocal motions with broken spatial symmetries 
for which the net translational velocity is proportional to Ref, for a certain a > 0, 
thus demonstrating a continuous breakdown of the scallop theorem for any Rcf > 
(see also [l3])- In this paper, we study the breakdown of the scallop theorem with 
Rcp in the presence of spatial asymmetries. While Lauga j3] studied the limit of 
{Re^, Rcp} <C Rcf <C 1, here we consider dense swimmers (pp/p 3> 1 and Pp/p ^ A/a) 
for which {Re^,Ref} -C Rcp <^ 1. We consider therefore the limit of negligible fluid 
inertia and show that locomotion is possible nonetheless. As in previous studies of 
the dynamics of simple swimmers [e.g., 0, H, Q, the examples presented here do 
not necessarily resemble real organisms, but they are intended to illustrate the theory 
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and provide physical insight. Analogous to [2I, our examples show a breakdown of the 
scallop theorem for any arbitrarily small Rcp > 0, since the net translational velocity 
(or rotation rate) for all our examples is proportional to Rcp. 

In section [21 we derive a general framework, using Lorentz's reciprocal theorem 



id llTj , to describe the motion of a dense swimmer. In section |3l we apply our general 
framework to study four examples of spatially asymmetric, dense reciprocal swimmers 
and show that they experience directed motion for arbitrarily small values of Rcp. The 
nonlinear interaction mechanism responsible for locomotion in the four examples is 
summarized in section |H 

2. Reciprocal theorem for a dense swimmer 

Consider the limit of a dense swimmer, {Rei^,Re{} <C Rcp -C 1, introduced in section 



[H If the fluid inertia is neglected, Lorentz's reciprocal theorem [16|] can be applied to 



relate the instantaneous dynamics of the swimmer to those of a towed, rigid body of the 



same shape 17|. Let the rigid body be towed with a constant force, F, plus a constant 
torque about a point P in the body, Lp. Assume the swimmer and the rigid body to be 
immersed in fluids of the same kinematic viscosity. The surface velocities and stresses 
on the swimmer, {u, cr), are related to those on the rigid body, {u, &), by 



n ■ a- ■ u dS = / n ■ a ■ u dS, (1) 

S{t) Js{t) 

where S{t) is the instantaneous location of the swimmer surface and n is the normal 
vector to it. The surface velocity of the swimmer can be written as u = Up + ft x 
{r — r-p) + u' , where Up is the translational velocity of the reference point P, 17 is 
the angular velocity, r is the position vector with an origin fixed at the inertial frame, 
and u' is the deformational component of the surface velocity. Note that the values of 
u' depend on the choice of P. All time derivatives are referred to the inertial frame. 
The surface velocity of the towed, rigid body is it = TJp + O x [r — rp). For a dense 
swimmer of constant mass, m, and homogeneous density, balances of forces and torques 
require 

n ■ (7 dS = m (2a) 

5(i) 

^^^^ (r - rp) xin-a) dS = (r^ - rp) x (rn^^ + (26) 

where and are the position vector and velocity of the swimmer's centre of mass, 
G, and /q is the inertia tensor referred to G. Introducing fl2aj) and fl26|) into ([1]) results 
in 

-(Up + ftx (re - rp)) + f2 - F Up - Lp f2 

n a u' dS. (3) 

Sit) 
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This derivation is valid in the absence of a body force. In the presence of a homogeneous 
gravitational acceleration, g, the terms mdua/dt in fl2aj) . fl2fej) . and ([3]) should be 
replaced by m (dwc/dt — g)- In ([3]), the rigid body velocities, L/p and f2, are arbitrary. 
By alternately taking = and £/p = 0, ([3]) yields two differential equations for the 
translational and rotational velocities of the swimmer. Up and O. 

As an example, consider the application of ([3]) to a squirming sphere of radius a, 
which undergoes purely tangential deformation. Take P = G, which is assumed to 
remain at the sphere centre. The differential equation for L^g is obtained by choosing 
= 0, for which Lp = 0. Then, F = —Q-k^oUq, and n ■ & U = — 3/i/(2a), which 
yields 

where 5*0 represents the sphere surface. The differential equation for O is obtained by 
choosing L/q = 0, for which F = 0. Then, Lq = —Sir^a^O and n ■ a ■ {O x n) = — 3/i, 
which yields 

3 f , 

; \- Si = / n X u db. (5) 

15/i dt 8na^ Jg, ^ ^ 

Note that (jll) and ([5]) reduce to equations 4 and 6 of [17] in the limit of Rcp = 0. Note 
also that, for a reciprocal squirming deformation, the time averages of (jl]) and ([5]) yield 
zero net translation and rotation. The same conclusion is obtained from ([3]) for any 
reciprocal squirmer, since S(t) = So is constant, and F, Lp, and & remain constant 
over the period of motion for constant tip and O. Therefore, a dense reciprocal swimmer 
needs to undergo deformation normal to its surface in order to propel itself. 



3. Examples of dense reciprocal swimmers 

The general framework derived in the previous section is here applied to four specific 
swimmers: two unequal spheres moving along a straight line, two unequal spheres 
moving along a circumference, a scallop- like swimmer, and a deforming sphere. These 
simple examples illustrate reciprocal locomotion in the absence of fluid inertia; the goal 
is therefore to compute the net translational velocity of each swimmer as a function 
of Rcp. For simplicity, gravity is not accounted for in these examples; the inclusion of 
gravity is straightforward, and its effect is discussed at the end of section |H 



3.1. Translation of two unequal spheres 

Consider the swimmer shown in figure [H which consists of two unequal spheres. The 
spheres have radii a and b, which are different but of comparable magnitude, that is, 
P = b/a = 0(1). Note /5 7^ 1 is required to break the spatial symmetry and yield 
net motion. The time-dependent distance between the sphere centres, L{t), is large 
compared with the sphere radii, so that L = 0(a/e), where e = a/{L) -C 1, and the 
angular brackets () denote a time average over the period of motion. The two spheres 
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are able to exert equal and opposite forces on each other, so that L(t) is a prescribed 
periodic function of time with radian frequency uj. The spheres move with velocities 
W^{t) = W{t) + w'{t) and Wb = W{t) - w'{t), where w' = 1/2 dL/dt, and W{t) is 
the unknown translational velocity of the reference point, P, which is chosen at the 
midpoint between the sphere centres. The parameter of interest is the net translational 
velocity, {W{t)). 



A 




Figure 1. Geometry of the two unequal spheres moving along their line of centres 
(z-direction). 



The reciprocal theorem ([3]) relates the dynamics of interest to those of a system of 
two spheres separated a distance L being towed at a constant speed W along their line 
of centres. By computing the effect of each sphere on the other as that of an equivalent 
point force, the hydrodynamic resistance forces in 2;-direction on the rigid, towed spheres 



are 



16 



l_^ + 0(62) 



2L 



(6a) 
(66) 



Application of ([3]) to this swimmer yields 



w 



(7) 



Next, we non-dimensionalize (I7j) using a and 1/u; as the length and time scales. The 
non-dimensional variables and parameters are r = cut, (3 = h/ a, W = eW/{auj), 
and A(r) = eL/a. We expand W in powers of Re^ = p^a^uj/ii <^ e <^ 1, i.e.. 
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W = Wo + RcpWi + .... With this, ([7]) becomes 

i?epA(l + P') + . . . j + ^[(1 + /5)A - 3/36 + Oie')]{Wo + RW, + . . .) 

= (-1 + n + + P + Oie^)] A^. (8) 

The solution to 0(i?ep) is 

f.. l (-l + /3)AdA/dr 

and (VFq) = 0. Note that the average translational velocity would remain zero even 
if higher-order terms were retained, i.e., there is no net translation of the two-sphere 
system to 0(i?ep). In contrast, the solution to O^RCp) includes terms of the form 
A (d^A/dr^), which yield a non-zero average. For example, suppose the distance between 
the spheres varies as A = 1 + eAsinr, where A = A/ a = 0(1) is the non-dimensional 
amplitude. In this case. 



3/j(/3- l)(-/3^ + 3/3-1) , 
6(l + /?)2 



(Wi) = : , , " '-A' + O(e^). (10) 



Then, in dimensional form, the net translational velocity is 

{W) = Re,au;e'^^^^^^^^^^^{A/ar + 0{Re,e\Rel). (11) 

Note that (W) tends to zero as L ^ oo, reflecting the fact that the net translational 
velocity is due to the interaction between the spheres. The smaller sphere advances 
in front {{W) > 0) for 1 < /? < (3 + \/5)/2, and the larger sphere advances in front 
{{W) < 0) for /3 > (3 + V5)/2. 

3.2. Circular motion of two unequal spheres 

Consider the swimmer shown in figure [2], which is a two-sphere version of the rotator 



devised by Dreyfus et al [15|. As was the case in the previous example, the spheres have 
radii a and b, which are different but of comparable magnitude, that is, P = b/a = 0(1). 
Again, /3 7^ 1 is required to break the spatial symmetry and yield net motion. The 
spheres are constrained to move along the circumference of centre P (which is held 
fixed) and radius R, which is much larger than a, R = R/a ^ 1. The angle between 
the spheres, 26{t), varies periodically with radian frequency in a prescribed manner. 
Specifically, 9{t) = Oq + e9i(t), with e = a/{L), L is the distance between the sphere 
centres, and 0(6*0) = 0(6'i) = 0(1). As a result of the interaction between the spheres, 
there is an unknown rigid-body rotation with angular velocity f2(t), whose time average 
we want to determine. The spheres move with linear velocities Ug,{t) and f/b(^) and 
experience hydrodynamic resistance forces F}^{t) and F^{t) in the azimuthal direction. 

If both spheres were towed along the circumference at a constant angular velocity 
Cl, the azimuthal hydrodynamic forces on the spheres, taken positive in counterclockwise 
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Figure 2. The inertial rotator. The spheres are constraint to move along the 
circumference of centre P and radius R. 



direction, would be [16 



-1 H sme + O(e^) 

AR sin 9 8R ^ \ 

3a cos^ 6 3a „ ^, -i. 

1 H sin6' + 0(e'^) 

AR sm9 8R ^ ^ ; 



Introducing the expression of the prescribed motion into ( 112 and ( 11261) yields 



F^ = GnfiaRn 
F^ = QnfibRCl 



-1 + ^ (70 - 7ie^i 
-1 + 1^(70-716^1 + 0(6^)) 



(12a) 
(126) 

(13a) 
(136) 



where 70 = (3cos^6'o — 1)/ sin^^o and 71 = cos 6*0(2 + 3sin^6^o)/ sin^6'o. With fll3ap and 
(113 6p . application of ([3]) results in 



(a^ + b')R- 



An 
dt 



b^R 



d^- 

dt2 



PP3- 



-(a 
5^ 



RQ 

+ 69 



(a3 + 63)^ 



^Q-{F^R + F^K)Q 



{Ft - Ft)Rf^. 



(14) 



Next, we non-dimensionalize (JT^ by defining Rcp = p^a^uj / ii ^ e ^ 1, r = ujt, f] = b/a, 
R = Re/ a, and fl = fl/{eu) and expand f2 = fio + RepCli + . . .. Thus, 



9 



'-Rer^R 



+ 1) + 



/?9 + l 
iP' + 1) 



sin^ 00 + 0(e) 



dfic 



+ . . . 
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(15) 



The solution to 0(i?ep) is 




~ _ (/3 - 1) de, r 3 (3 70 
° {(3 + 1) dr [ ^ A{l + f3)R' 



6'id^6'i/dr^, which yield a net angular velocity. If we consider the case 9i = A/ (2i?) sinr, 
where A = A/ a is the non-dimensional amplitude of motion, the time-averaged non- 
dimensional angular velocity becomes 




It is noted that this non-dimensional velocity becomes zero when 6q = 7r/2, as this 
corresponds to the fore-aft symmetric configuration in which the spheres are diametrally 
opposed. The velocity becomes unboundedly large if R remains constant and ^ 0. 
In this limit, however, the spheres become close and f|T71) is invalid. According to f|T7|) . 
the smaller sphere advances in front for 1 < (3 < 1.74, while the larger sphere advances 
in front for /? > (3 + V5)/2 ^ 2.62. For 1.74 < (3 < 2.62, the direction of rotation 
depends on the value of Oq. 

In dimensional form, the net angular velocity is 



Note that the rotator reduces to the previous example of two unequal spheres moving 



sin 6*0 ^ 0, 7i ^ 8-R^, and the net linear velocities predicted by (ITT]) and (fT8l) agree. 

3.3. Translation of a scallop-like swimmer 

We study the scallop-like swimmer schematized in figure [31 The swimmer consists of 
two circular cylinders of radius a and length L ^ a connected by a hinge. The angle 
between the two cylinders, 26{t), is a prescribed periodic function of time with radian 
frequency u. The swimmer undergoes an unknown translational velocity, W{t), whose 
time average we want to determine. 

By neglecting hydrodynamic interactions between the cylinders, the hydrodynamic 
forces exerted on each cylinder are related to the velocities of the cylinder by the 
resistance coefficients. The resistance coefficients for translation perpendicular and 



(Q) = ujeRep{Qi) + 0{e^Rep, i?ej). 



(18) 



along their line of centres for R 



oo with constant e = l/(2i?sin6'o). In this limit. 
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Figure 3. Scallop-like swimmer consisting of two circular cylinders of radius a and 
length L connected by a hinge at P. 



parallel to the cylinder axis are R 
of ([3]) in the z-direction yields 

Lde 



2i?|| = An^L/\n{L/a) [e.g., [l8|. Application 



m 



dW d_ 
dt 



sin 9 



W 



2 dt 

+ (-R±W sine ]W sine 



-RiiW cose] W cose 



-R±W sine 



Lde 
2"dt' 



(19) 



where m = p^TTo?L is the mass of each cylinder. This equation is made non-dimensional 
by using a and l/cj as length and time scales and by defining a non-dimensional time 
r = ut, cylinder length \ = L/a, and velocity W = W/ {uL). Introducing these scalings 
into (|T9|) and expanding W in powers of Re^ = p^a^uj/ ii ^ e ^ 1 result in 



RCr 



dWo 



+ ... + 



In A 



(l + sin^^) (Wo + Re^Wi + ..^ 



Rep d 
2~d7 

The solution to O(-Rep) is 



2 . Ae 

sm e—. 

dr 



sin e de 



(20) 



(21) 



1 + sin^ e dr 

and (Wo) = 0. Let us assume small oscillations of the form 6'(r) = 6'o+6'i(r) = 6'o+e sin r 
where e <^ 1. Then, the solution to O(-Rep) has a non-zero average given by 

1 



(Wi) = -e^ln(A)/(^o) + 0(e^ 



where 



f{0o 



sin? 6*0 cos'^ 6*0 



(22) 



(23) 



(l + sin^ ^o)'' 

The dimensional net velocity is {W) = RepUjL{Wi) + O(-Rep). The function /(6'o) is 
represented in figure HI Since /(6'o) > for all eo e (0, 7r/2), {Wi) > 0, and the 
swimmer moves with the hinge advancing in front. The net velocity is maximum for 
^0 = arccos([-7/2 + y73/2]i/2) ~ 28.5°. 
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O.lr 




Jl/8 Jl/4 3JI/8 Jl 



^0 

Figure 4. Dependency of the net translational velocity on the average angle of opening 
of the scallop- like swimmer, ^o- 

3.4- Translation of a deforming sphere 

A quasi-spherical swimmer of radius a executes a small, periodic deformation, so that 
its radius is given by R(t,9) = a(l + e^ao(^)-Po(cos ^) + eai()f:)Pi(cos6')), as shown in 
figure [5l Here, Pq = 1 and Pi = cos 6 are the Legendre polynomials of orders and 
1, respectively, and e ^ 1. The coefficient of Pq is chosen as ao = —ai/3 + O(e^), 
so that the particle volume remains constant and equal to Ana^/S. The deformation 
is axisymmetric with respect to the z-axis, and it is such that the particle density, pp, 
remains homogeneous and constant. 




Figure 5. Undeformed (dashed line) and deformed (solid line) geometry of the 
swimmer. The deformation is axisymmetric with respect to the z-axis. Not to scale. 

To compute the translation of the sphere, we apply the 2;-component of ([2]) to this 
swimmer, defining the reference point P as the centre of the undeformed sphere. To 
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evaluate the first term on the left-hand side of ([3]), we write the velocity of the centre of 
mass, G, as uq = up +Mq, where u'q is the velocity of G relative to P. The acceleration 
of G relative to P is 



dt 



d^ai 



(24) 



directed along the z-axis. The normal vector to the deformed sphere's surface in 
spherical coordinates (r, 6, 0) is 



n 



+ (di?/d6»)2/_R2]l/2 

{dR/de)/R 



\ 



[i+{dR/de)2/-R2] 




T72 



f 1 - ^a\Vll2 \ 
taiVi - e^alPiVi 




(25) 



where we have used Lighthill's 19| definition, 



2VT^dP„,(r/) 



(26) 



n{n + 1) drj 

with 7] = cos 6. To evaluate the two surface integrals in ([3]), we must calculate the flow 
around the rigid, deformed sphere past a uniform flow, U. Any axisymmetric flow can 
be described by the streamfunction 



n=0 

where An, Bn, C„, and Dn are constants, and d(5„(?7)/d?7 = Pn{f]) [^-g-, 
and azimuthal velocity components, Ur and ue, are given by 

1 dijj 

1 di) 

Ug = - 



20 



Ur 



(27) 
The radial 

(28 a) 
(285) 



r^/l - r/2 dr 

In our problem, in which the reference system is attached to the swimmer, Ur and 
ug must tend to the ambient velocity, of magnitude ?7, as r — oo. In addition, the 
solutions with n = are unadmissible, because the corresponding azimuthal velocities 
are undefined at = 1 (9 = 0). With these constraints, the general solution of the flow 
field is 



Ur 



Ue 



n=l 



rn+2 



t/sin^ + 5: + 

n=l ^ 



Pn 



as obtained in previous studies of spherical swimmers pj|, 121 
determined using the boundary conditions, = = 0, at r 
eaiPi). These constants can be expanded in powers of e: 



(29 a) 

(296) 

The constants are 
R = a{l + e^ttoPo + 



(30 a) 



fe=0 
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(306) 

(31a) 
(316) 
(31c) 
{31 d) 

while all other constants are at most of O(e^). The value of the pressure on the deformed 
sphere's surface is 





= Y.''B^\ 


is 






A, 


= U 




Bi 


= u 


' 1 2/^3 


A, 


= U 




Bs 


= U 





10 



piR) 



9 



9 



4 



+0(e^)] . 

The radial and tangential stresses at r = i? are 



a. 



rr\r=R 



fiU 

a 



--Pi - e«i3P/ 



33 



27 



( -Pt + ( 7T«o - -at ] Pi + ^«f^3 ) + 0(e 



237 



^re\r=R 



/if/ 

a 



-Vi - e«i6PiVi 

(|«?PiVi + (-^ao - ) V^i 
+0(e^ 



153 

To" 



am 



<^ee\r=R 



a 



--Pi + eai6Pi2 + 0(e2) 



Finally, the two surface integrals in ([3]) are 

/ cos 9(Trr — siu OarO ^ 



n ■ & ■ z dS 



S{t) 



n 



cos 6are — sin dcxs 
cos ^cr^^ — sin 6as 

7 



(27rp2 sin eae) 



J 



1 + eUao- —af + 0(e 



/ 

Js{t) 





/ 


CXrr 


\ 




& ■ u' dS = n ■ 

Jo 




CXre 




a ^ 




\ 


0'r4) 







15 



2 daQ 



dt 



+ ePi^^ {27iR^smed9) 



-2e^ + 
dt 



2d«i dai dao 
'2ai—^ 2ao— 4q;i 



dt 



dt 



dt 



(32) 



(33 a) 



(336) 

(33 c) 



(34) 



0(e^ 



.(35) 
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Introducing these results into yields, in non-dimensional form, 

U 



4„ dU 



i + f=(-5"i-^"?)+0(«') 

^dai 2/^2 ^ \ dai 16^ d^ai ^, o, 
= + ' (-2«? + 2«i) ^ - Y^R^P^ + 0(e ). (36) 



where r = cut, f/ = euaU. Next, we expand U = Uq + Re^Ui + . . . and assume 
i?ep -C e ^ 1. The solution to O(-Rep) is 

1 dai , / 22 9 2 \ dai ^ , , 

For any periodic ai(r), the time average of this expression is zero. Moreover, the time 
average is zero to all orders of e. 

The solution to 0{Re^) includes terms of the form e^(dai/dr)^ and e^Q;i(dai/dr)^ 
which, for a suitable choice of ai(t), yield a non-zero average. For instance, for 
a\ = sin^(a;t) = sin^ r, the time-averaged non-dimensional velocity is 

{U,) = -^e' + 0{e% (38) 
or, in dimensional form, 

{U) = -^Re^e^iua + 0(e^i?ep, i?ej). (39) 
Note that the minus sign indicates a net translation in the (— 2;)-direction. 

4. Conclusion 

We have presented several examples of dense reciprocal swimmers that are able to propel 
themselves, even in the absence of fluid inertia, at arbitrarily small values of Rcp. The 
existence of a net swimming velocity arises from a nonlinear interaction between the 
oscillatory particle inertia and the oscillatory drag in the presence of spatial asymmetries. 
The existence of this nonlinear interaction can be inferred from the reciprocal theorem, 
(IH]). In the absence of rigid-body rotation, (l3l) reads 

(^^) ■ - ^ ■ = / ri & u'dS. (40) 

F, the hydrodynamic drag on the towed, deformed body, scales as {fj,aUp){l + A/aa(t)). 
Here, (fiallp) is the magnitude of the hydrodynamic drag on the towed, undeformed 
body, A is the amplitude of oscillation of the swimmer surface, a{t) is the non- 
dimensional oscillation, and a ^ Ais the size of the undeformed swimmer. The integral 
on the right-hand side is typically of magnitude {fiaUp){Ki + A/a K2a(t))u' , where 
u' = A da(t)/dt is the magnitude of the deformational surface velocity u', and Ki and 
K2 are coefficients that account for the change of magnitude and direction of u' over 
the swimmer surface. Also, ~ f/p ~ Auj. Introducing these scalings into fHOj) yields 

R^p^ - (1 + ea{r)) Up ~ {K, + eK2a{T)) ^, (41) 
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where Up is the non-dimensional translational velocity, t = ut is the non-dimensional 
time, and e = a/A <^ 1. By expanding Up in powers of Rcp, the equation to O^Re^) 
yields 

f/p,o r^-[K^ + e {K2 - Ki) a(r) + 0(6^)] — , (42) 

and (Upfi) = 0, in agreement with the scallop theorem. Next, the equation to O(-Rep) 
yields 

Up,i ~ ^ [1 - ^«(^) + 0(e^)] . (43) 

Thus, in the presence of certain spatial asymmetries, the right-hand side of (H3l) contains 
a term of the form 



d^a , , /da^ ^ 



{K2 - 2K,) a(r)— + (7^2 - i^i, , 



(44) 



For an appropriate choice of the reciprocal oscillation ^(t) (for example, a sinusoidal 
function), this term gives raise to a non-zero net velocity to O(-Rep). Accordingly, in the 
four examples presented in this paper, the net translational velocity is proportional to 
Rcp. As discussed in section |2l a necessary condition for the existence of this non-zero 
velocity is that the dense swimmer undergoes deformation normal to its surface, while 
pure squirming does not result in locomotion. 

In the previous analysis we have not included the effect of gravity, which needs to 
be considered in practical applications. For A ~ a and e = 0(1), we expect the net 
translational velocity due to the reciprocal motion to scale as 

(Up) ~ auRcp, (45) 

while a typical sedimentation speed is 

U, ~ (46) 



Thus, the dense swimmer needs to oscillate at a frequency uj > 0{^Jg/a) in order to 
overcome gravity. 

In conclusion, our examples show the existence of a net translation or rotation at any 
arbitrarily small Re^ > for certain spatially asymmetric, dense reciprocal swimmers. 
This demonstrates the breakdown of the scallop theorem for a case in which the flow 
around the swimmer is governed by the Stokes equations, that is, in the absence of fluid 
inertia. While we have studied the limit of dense swimmers, for which pp ^ p, in many 
biological applications Pp ^ p and the unsteady fluid inertia is expected to be comparable 
to the particle inertia. The extension of this study to the limit of Re^ <^ Re^j ~ Re^ 
will be the subject of future work. 
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